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ABSTRACT 

We have re-analyzed continuum and recombination lines radio data available in the literature in order to 
derive the luminosity function (LF) of Galactic HII regions. The study is performed by considering the first 
and fourth Galactic quadrants independently. We estimate the completeness level of the sample in the fourth 
quadrant at 5 Jy, and the one in the first quadrant at 2 Jy. We show that the two samples (fourth or first quadrant) 
include, as well as giant and super-giant HII regions, a significant number of sub-giant sources. The LF is 
obtained, in each Galactic quadrant, with a generalized Schmidt's estimator using an effective volume derived 
from the observed spatial distribution of the considered HII regions. The re-analysis also takes advantage of 
recently published ancillary absorption data allowing to solve the distance ambiguity for several objects. A 
single power-law fit to the LFs retrieves a slope equal to -2.23±0.07 (fourth quadrant) and to -1.85±0.1 1 (first 
quadrant). We also find marginal evidence of a luminosity break at L^nee = 10 23 45 erg s -1 Hz" 1 for the LF in 
the fourth quadrant. We convert radio luminosities into equivalent Ha and Lyman continuum luminosities to 
facilitate comparisons with extra-galactic studies. We obtain an average total HII regions Lyman continuum 
luminosity of 0.89 ±0.23 x 10 53 sec _1 , corresponding to 30% of the total ionizing luminosity of the Galaxy. 
Subject headings: luminosity function, HII regions, radio continuum 



1. INTRODUCTION 

Extensive studies have investigated the LF of HII regions in 
external galaxies. These, mostly based on Ha observations, 
have established that: 1) the slope of the LF is correlated 
with Hubble type, with early-type galaxies having a steeper 
slope than late-type (Kennicutt, Edgar & Hodge 1989, here- 
after KEH89); 2) the LF becomes significantly flatter at low 
luminosities (KEH89; Rand 1992; Walterbos & Braun 1992; 
Rozas et al. 1996, 1999, 2000; Bradley et al. 2006). Such ob- 
servational facts have been interpreted in terms of evolution- 
ary effects (Oey & Clarke 1998). Variations of the LF in spiral 
arm versus inter-arm regions are controversial: some studies 
claim that the LF of interarm HII regions is steeper (KEH89; 
Banfi et al. 1993; Rand 1992), but others do not support this 
conclusion (Rozas et al. 1996; Knapen et al. 1993, Knapen 
1998). 

Remarkably, there is a surprising paucity of similar stud- 
ies regarding the HII regions population of our own Galaxy. 
This is partly due to the fact that HII regions lie on the Plane, 
where Ha is heavily obscured by dust, therefore limiting the 
observation to HII regions located within a radius of a ~ kpc. 
Radio surveys could in principle circumvent this problem but, 
in the past, only a few have uniformly targeted the Galactic 
Plane region with sufficient angular resolution and sensitiv- 
ity to resolve individual sources. Likely the near future will 
see a breakthrough in this direction, with the completion of 
extensive radio surveys carried out through interferometric 
techniques, such as the International Galactic Plane Survey 
(IGPS) at 1.4GHz and 408 MHz, or MAGPIS (White et al. 
2005), at 5 and 1.4 GHz. Kerton et al. (2007) have pub- 
lished a catalog of Galactic extended sources from the Cana- 
dian Galactic Plane survey (Taylor et al. 2003), part of IGPS. 
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However, despite a remarkable work of detection and extrac- 
tion, the sources have not been classified yet. 

The reference works for the LF of Galactic HII regions are 
Smith & Kennicutt (1989) and McKee & Williams (1997). 
They both rely on the list of HII regions compiled by Smith, 
Biermann & Mezger (1978, hereafter SBM). This list is a 
compilation of 93 HII regions and it is characterized by im- 
portant caveats. For example: the Southern sources in the 
list (49 out of 93) are mainly from the 5 GHz recombination 
lines survey by Wilson et al. (1970). These authors claim that 
their survey is complete for sources having a peak antenna 
temperature of 1.3 K, or greater, at 5 GHz. However, this 
statement is based on the continuum data which have been 
obtained through a variety surveys, carried out with different 
frequencies and angular resolutions. For instance, the Thomas 
& Day (1969a,b) 2.7-GHz survey has a spatial resolution of 
8.2', while the Hill (1968) 1.4-GHz survey is characterized by 
a resolving power of 14'. Most importantly, the completeness 
levels of these continuum surveys appears to be unknown. 

Therefore, Wilson et al. cannot be certain that all sources 
above the 1-K thereshold have been picked up, in the first 
place, by those surveys. This means that the recombina- 
tion lines survey has an intrinsic bias, with important con- 
sequences on the completeness for giant and super-giant HII 
regions, as we will discuss in § 5. 

Complementary studies have focused on the LF of ultra- 
compact HII regions (Comeron & Torra 1996; Casassus et al. 
2000). The work by Comeron & Torra is probably the only at- 
tempt found in the literature to take into account the complex 
spiral-arm geometry of our Galaxy. However, since it relies 
only on the IRAS point-source catalog (PSC) and lacks any 
kinematic information on the sources, it purely consists in a 
multi-parameter fit analysis, which has obvious limitations. 

For the reasons described above, we have decided to look in 
the literature and search for alternative data to the SBM list in 
order to obtain an independent derivation of the LF of Galactic 
HII regions. We are in fact convinced that archival data are 
widely and unjustifiably under-used in this field and, until new 
catalogs of HII regions are made available, for instance, by the 
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MAGPIS or the IGPS consortia, these archival data should be 
fully exploited. This is even more true given the new release 
of ancillary data, such as of those on HI absorption, which 
make it possibile, for many of these sources, to break, for the 
first time, the distance degeneracy. 

In the following, we will consider two definitions of the LF 
denoted, respectively, by n(L)dL and (f>(L)dlogL. The for- 
mer gives the number of sources per unit volume with lumi- 
nosity in the range [L, L + dL] and is usually parameterized 
as n(L)dL oc L~ a dL; the latter gives the number of sources 
per unit volume with luminosity in the logarithmic interval 
[logL, logL+fflogL] and is represented as <p(L) dlogL oc 
L~ a+l dlogL. 

2. SAMPLE SELECTION 

The vast majority of catalogued HII regions is located in 
the first (0° < I < 90°) and fourth (270° < / < 360°) Galactic 
quadrants. 

The most recent, extensive radio recombination lines sur- 
veys in the Galactic Plane are the Caswell & Haynes (1987, 
hereafter CH87) survey in the fourth quadrant, and the Lock- 
man (1989) survey in the first quadrant. 

The CH87 survey has covered the region 210° < I < 360°, 
\b\ < 2°, in the fourth quadrant. With an angular resolution 
of ~ 4.1', it represents the follow-up of the continuum 5 GHz 
survey by Haynes et al. (1978, 1979), which has been car- 
ried out at the same angular resolution as the recombination 
lines survey. The CH87 survey consists of HI 09a and HI 10a 
recombination lines measurements of most of the sources (as 
stated by the same authors) with peak brightness temperature 
exceeding 1 K at 5 GHz, equivalent to a completeness limit, 
in peak flux density, of 1.3 Jy. CH87 provide, for each HII 
region in their surveyed region, the measured recombination 
line velocity and electron temperature, as well as the angu- 
lar diameter and the integrated flux density as observed in the 
continuum 5 -GHz survey. We point out that these flux den- 
sities and diameters have been derived by Haynes et al. by 
measuring by hand contours of surface brightness. The same 
method has been applied to the sources in the SBM list lo- 
cated in the Southern hemisphere, as described in Wilson et 
al. (1970). 

The Lockman (1989) survey, with a spatial resolution of 
3', is a follow-up of the 4.8 GHz Altenhoff et al. (1979) 
continuum survey, made with comparable angular resolution 
(2.6')- It consists of observations of the H87a, H88a and 
H85a hydrogen recombination lines, and it is complete down 
to 1 Jy, in peak flux density, for 0° < / < 60°, while for 
I > 60°, the completeness of the survey is less uniform. In- 
terestingly, the Lockman survey has covered, with equal an- 
gular resolution, the same region of the sky of the HI 10a and 
formaldehyde (H 2 CO) absorption line survey by Downes et 
al. (1980, hereafter DWBW). DWBW also claim to be com- 
plete to a peak flux density of 1 Jy, although Lockman re- 
ports that <~ 200 sources are observed for the first time. This 
fact could be explained with misidentification problems in the 
DWBW survey. For the Lockman sources, flux densities and 
angular diameters are taken from the 4.8 GHz Altenhoff et 
al. (1979) survey, and these have been measured through au- 
tomated gaussian fits. The same technique has been used to 
derive flux and angular sizes for the sources in the SBM com- 
pilations located in the Northern hemisphere, as discussed in 
Reifenstein et al. (1970). 

Instrumental systematics can potentially degrade in a sig- 
nificant way continuum measurements, especially with the 



technology available 20 or 30 years ago. One of the main 
problems is represented by 1 // noise which can cause fluctu- 
ations of the baselines and make the measured flux densities 
unreliable. The solution is typically to scan the sky at a fast 
rate. The scan rates of the Haynes et al.'s and Lockman's sur- 
veys are 2.5°/min and 80'/min, respectively. In comparison, 
the sources above 1 K at 5 GHz and selected by Wilson et al. 
for their recombination lines survey, have been re-observed in 
the continuum by Goss & Shaver (1970) at a rate of l°/min, 
meaning that the uncertainties in the flux densities quoted by 
CH87 and Altenhoff et al. (1979), due to effects such as 1 // 
noise, are lower or at least of the same order as the values 
reported by SBM. 

We note that Paladini et al. (2003; hereafter Paper I), com- 
piled a radio catalog of 1442 Galactic HII regions for which 
angular diameters and flux densities at a reference frequency, 
namely 2.7 GHz, are given. This is the so-called Master Cat- 
alog, in which, for each source, the quoted values at the refer- 
ence frequency are obtained by taking the weighted mean over 
all the available data reported in the literature. This proce- 
dure also requires to estimate measurement errors when these 
are not provided by the authors of the original surveys. For 
~ 800 of the HII regions in the catalog, radio recombination 
lines have been measured and the mean line velocities, line 
and electron temperatures are reported in the Master Catalog. 
Such a data set, combined with the Fich, Blitz & Stark (1989, 
hereafter FBS89) rotation curve, has been used in Paladini, 
Davies & De Zotti (2004; hereafter Paper II) to compute, for 
a subset of 575 sources 5 , galactocentric and solar distances. 

We could have in principle used the Paladini et al.'s cat- 
alog to perform our analysis of the LF. However, as men- 
tioned above, the Master Catalog contains weighted averages 
of measurements obtained with different instruments, at vari- 
ous frequencies and angular resolutions. This fact per se can 
introduce complications in, for example, assessing the actual 
completeness of the sample. 

On the contrary, the CH87 and Lockman's data bases 
present the advantage (over the Paladini et al.'s catalog but 
also over the SBM list) of being uniform samples, obtained by 
observations of large areas of the sky with a single instrument, 
characterized by a given sensitivity and angular resolution. 

3. COMPLETENESS OF THE SAMPLES 

CH87 observed 316 HII regions in total. The authors state 
in the paper that they have likely observed most of the sources 
with S p > 1 .3 Jy. The first step is then to remove from the sam- 
ple the sources below this threshold, i.e. since CH87 quote, 
for each source, the integrated flux at 5 GHz, 5,, this corre- 
sponds to exclude the sources for which: 
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where 6>hii is the angular diameter of the source, and Ob is 
the angular resolution of the survey (4.4'). The application 
of the above criterium leaves us with 238 HII regions. At 
this point, we still need to establish the actual completeness 
of the survey in terms of flux density. For this purpose, we 
consider the observed differential counts (see Fig. 1) of the 
sample and we fit to this, along the lines of Paper II, a simple 

5 Only the sources with a measured velocity |V/„-| > 10 km/s are consid- 
ered in the analysis. 
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function of the type S' 1 , representing the scaling of the num- 
ber of sources, as seen from an observer located at the Sun, 
if these were distributed along a disk. This approximation is 
valid within the range 30 Jy > Sschz > 5 Jy. For Sschz > 
30 Jy, the distribution steepens and is better represented by 
a function S~ 3 / 2 , appropriate for a uniform three-dimensional 
distribution, roughly corresponding to nearby sources. Errors 
have been computed, in each flux bin, by applying the formu- 
lae by Gehrels (1986). Specifically, the upper limits are given 
by: 

(2) 
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while the lower limits are: 
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where S = 1, j3 = and n is the number of sources per bin. 
By removing from the sample the objects with an integrated 
flux < 5 Jy, we are left with 162 sources. In the following, 
we will derive the LF of the fourth quadrant also for flux cut- 
offs different (and above) the cut-off value obtained from the 
the analysis of the differential counts. These additional LFs 
will be used for comparison with the one obtained for the ref- 
erence 5-Jy cut-off, in order to estimate the potential error 
introduced by the assesment of the completeness level for the 
CH87 sample. 




FIG. 1. — Differential source counts for the CH87 sample. The fits for a 
disk (dashed-dotted line) and spherical (dashed) distributions are shown. The 
completeness level is denoted by the dashed vertical line. 

The Lockman (1989) sample has 256 sources within longi- 
tude 2° < I < 60°. After removing the sources with peak flux 
density S p < 1 Jy, or for which: 
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(taking h = 2.6'), we are left with 1 15 sources. In this case, 
the source counts are well fitted by a model <~ S~ l in the flux 
range 30 Jy > Ssghz — 2 Jy, while for Ssch z > 30 Jy a function 
5~ 3//2 is a better representation of the data. The exclusion of 
sources below the 2-Jy thereshold leaves 93 objects. Along 
the same lines of the CH87 sample, we will also compute the 
LF in the first quadrant for completeness cut-offs above the 
one estimated from the differential counts. 

4. EXCLUSION OF SOURCES WITH \V LSR \ < 10 KM/S 

Sources with radial velocities | Vlsr| > 10 km/s can, in prin- 
ciple, be affected by non-circular streaming motions, but, as 
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FIG. 2. — Differential source counts for the Lockman (1989) sample. The 
fits for a disk (dashed-dotted line) and spherical (dashed) distributions are 
shown. The completeness level is denoted by the dashed vertical line. 

noted for instance by Quireza et al. (2006a), an accurate quan- 
tification of this effect is not possible given our current knowl- 
edge of the Galactic velocity field. For this reason, throughout 
the text, we have assumed that the circular motion hypothesis 
is generally valid. 

On the contrary, sources with measured radial velocities 
< 10 km/s are known to be significantly contaminated by pe- 
culiar velocity components, and distances derived by using 
these velocities can be highly unreliable. 

Therefore, from each sample (CH87 and Lockman (1989)), 
we remove the sources with |Vlsr| < 10 km/s. For the refer- 
ence cases of flux cut-offs of 5 and 2 Jy, this corresponds to 
remove 23 objects in the fourth quadrant and 5 objects in the 
first quadrant. 

In § 8.3, we will examine the effect, on the derived LF, due 
to the removal of the low- velocity sources. 

5. RADIAL AND SOLAR DISTANCES 

For each source in the CH87 and Lockman (1989) sam- 
ples, we compute the galactocentric and solar distances, R 
and D, respectively, by combining the measured hydrogen re- 
combination line velocity with the rotation curve by FBS89. 
While in Paper II we used Rq = 8.5 kpc, we have adopted here 
a somewhat smaller value, Rq = 8 kpc, favoured by recent de- 
terminations (Eisenhauer et al. 2003, 2005). FBS89 provide a 
useful tabulation of the rotation curve for several values of Rq 
(see Table 4 of their paper), including Ro = 8.0kpc. In § 8.2, 
we will show that the adoption of rotation curve models other 
than FBS89 has a very minor impact on the derived LF. 

As for the solar distances, part of the HII regions in both 
samples are affected by the well-known distance degeneracy 
problem. These are the sources within the solar circle for 
which two solutions are consistent with the observed radial 
velocity. In this case, we proceed as in Paper II, using auxil- 
iary (optical or absorption) data, when available, to break the 
degeneracy. 

For the fourth quadrant (CH87 sample), with respect to Pa- 
per II, we have complemented the list of references of auxil- 
iary data with the work on HI absorption by Fish et al. (2003), 
Nord et al. (2006) and Quireza et al. (2006b). 

When sources lack auxiliary data, we make use of the 
luminosity-physical diameter correlation found in Paper II 
and, in particular, of eq. (7) in that paper, which, for the 5 
GHz case takes the form: 
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with C = 1 .2767 x 10" 19 , a = 32.4 and b = 0.86. 6 For the ref- 
erence case of a flux cut-off at 5 Jy, we have applied eq. (5) to 
61 out of 139 sources, corresponding to 44% of the total. 

We emphasize that the relation above is used to choose be- 
tween the near and far kinematic solutions for a given mea- 
sured line velocity, not to compute a heliocentric distance 
from the observed angular diameter and flux. Following 
Paper II, we performed on the CH87 sample of sources a 
test to cross-check whether the application of eq. (2) to the 
sources with auxiliary data assigns the correct solution. The 
test shows that the luminosity-diameter correlation allows to 
choose the correct solution in ~ 80% of the cases. 

Fourth quadrant 
50 b J 




4 5 6 7 8 9 10 

R [kpc] 



FIG. 3. — Radial distribution of HII regions in the 4th Galactic quadrant. 
Sources are from CH87, and their flux density is above 5 Jy at 5 GHz. Bins 
of AS = 1 kpc. The solid line denotes the distribution obtained when sources 
with \Vlsr\ < 10 km/s are left in the sample, while the dashed line corre- 
sponds to the case when these sources are removed. 

In the first quadrant, 79 out of 88 of the Lockman sources 
above the reference flux cut-off at 2 Jy lie inside the solar cir- 
cle. The first step is to make use of auxiliary (optical and ab- 
sorption) data to solve the kinematic distance ambiguity. As 
for the CH87, we complement the list of references of Paper II 
with the HI absorption data by Fish et al. (2003), Kolpak et al. 
(2003), Quireza et al. (2006a), Anderson & Bania (2009), as 
well as with the formaldehyde (H2CO) absorption data pro- 
vided by Araya et al. (2002), Watson et al. (2003), Sewilo 
et al. (2004) and the same DWBW. However, as noted by 
Anderson & Bania, HI absorption methods are usually more 
robust than H2CO to break the distance degeneracy, given that 
HI is more ubiquitous than H2CO. For this reason, whenever 
a source has both H2CO and HI absorption data, we favor the 
solution associated to HI. Following these criteria, we are able 
to solve the distance ambiguity for 70 objects. To the remain- 
ing 9 (10% of the total) sources we have applied eq. (5). 

The radial distributions of the HII regions above the flux 
cut-offs of 5 and 2 Jy are illustrated in Fig. 3 and Fig. 4, 
for the fourth and first quadrant, respectively. Note that the 
figures show the radial distribution both when sources with 
\Vlsr\ < 10 km/s are removed or left in the samples. Clearly, 
the removal of these sources causes, in the fourth quadrant, an 
artificial dip in the distribution for R ~ 8 kpc. 

Remarkably, and regardless of the exclusion of the low- 
velocity sources, the radial distribution appears to be differ- 
ent in the two quadrants. As emphasized by several authors 
(e.g. Lockman 1981, 1989; CH87), this is a consequence of 

6 Note that, in Paladini et al. (2004), eq. (7) is missing the factor 1.2767 X 
10" 19 and the frequency v is in wrong units, i.e. GHz instead of Hz. 



the fact that HII regions, as massive stars in general, trace the 
spiral structure of the Galaxy and, as such, reflect in their spa- 
tial distribution the asymmetry of the arms geometry as seen 
by an observer located at the Sun position. 
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FIG. 4. — Radial distribution of HII regions in the 1st Galactic quadrant. 
Sources are from Lockman (1989), and their flux density is above 2 Jy at 4.8 
GHz. Bins of AR = 1 kpc. As for the fourth quadrant: the solid line denotes 
the distribution obtained by including all sources; the dashed line denotes the 
one corresponding to the removal of sources with \Vlsr\ < 10 km/s. 



6. COMPOSITION OF THE SAMPLES 

The SBM compilation used in Smith & Kennicutt (1989) 
and McKee & Williams (1997) is supposedly complete, over 
the whole Galaxy, for super-giant HII regions, and for our 
side of the Galaxy, for giant HII regions. SBM define super- 
giant HII regions those having a Lyman continuum luminosity 
N c > 50 x 10 49 s -1 , while giant HII regions are the ones for 
which 4x 10 49 s" 1 > N c > 50x 10 49 s _1 . These definitions 
are provided taking as a reference the luminosity of Orion 
for Rq = lOkpc and T e = 10000 K. Since in our derivation of 
the radial and solar distances we use Rq = 8. kpc and T e ~ 
7000 K, we rescale, following McKee & Williams, the Lyman 
luminosities by a factor 0.75 which takes into account both 
corrections. Accordingly, the new limits for super-giant and 
giant HII regions become, respectively, N c > 37.5 x 10 49 s" 1 
and3x 10 49 s" 1 >N C > 37.5 x 10 49 s~ f . In addition, sub-giant 
HII regions will be the sources with an ionizing luminositiy 
less than 3x 10 49 s" 1 . 

We use the flux densities tabulated by Haynes et al. (1978) 
for the CH87 sources and those by Altenhoff (1978) for the 
Lockman sample, combined with the solar distances assigned 
as described in the previous sections, to compute for each HII 
region in the fourth and first quadrant its ionizing luminosity, 
following the Condon's (1992) relation: 



Q = 6.3 x lO^photonss" 1 • (6) 

A GHz J 10 27 ergs" 1 Hz -1 

where, for the electron temperature, we use T e = 6200 K for 
the sources in the fourth quadrant, and T e = 7400 K for the 
objects in the first quadrant. These numbers represent mean 
values of the electron temperatures of the sources in the CH87 
and Lockman's samples. 
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TABLE 1 

Composition of CH87 and Lockman samples 



Type 


CH87 


Lockman (1989) 


sub-giant HII regions 


74 


51 


giant HII regions 


62 


35 


super-giant HII regions 


3 


2 



NOTE. — Quoted values are for the reference cases of a flux cut-off at 5 
and 2 Jy, for the CH87 and Lockman (1989) samples, respectively. 

Applying the definition of sub-giant, giant and super-giant 
HII regions given above, we obtain the composition of the 
CH87 and Lockman samples provided in Table 1 . 

Clearly, in both quadrants we have a significant contribu- 
tion from sub-giant HII regions. This is important given that, 
as noted by McKee & Williams (1997), most of the ioniz- 
ing luminosity of the Galaxy is likely contributed by small 
HII regions, rather than by the super-giant and giant popula- 
tions. Furthermore, we note that, for the fourth quadrant, the 
SBM list quotes 49 between giant and super-giant HII regions, 
while the CH87 sample includes 65 in total. Likewise, in the 
first quadrant, the SBM compilation has 33 giant and super- 
giant HII regions in the region 2° < I < 60°, compared to 37 
for the Lockman (1989) sample. This shows that the sam- 
ples considered here are likely more complete than the SBM 
catalog, as mentioned in the Introduction. 

7. THE LF FOR AN INHOMOGENEOUS SOURCE 
DISTRIBUTION 

One of the limitations of previous works which have inves- 
tigated the LF for Galactic HII regions is that they did not 
allow for a realistic geometric distribution of the sources. In 
this section we describe the approach adopted for our analy- 
sis. Under the assumption that the LF is independent of posi- 
tion, the number N(L) of HII regions with luminosity L within 
dlogL detected above a flux threshold S is related to the LF, 
4>(L)d log L, by: 



(7) 
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with R = (Rl + D 2 -IDRqcosI) 1 ! 2 . In the above expression 
n(R) and p(z) describe the radial and vertical density profiles, 
respectively, and the multiple integral on the right-hand side 
is the effective volume within which a source of luminosity L 
is included in the catalog, i.e.: 
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where £> max (S|L) = (L/4irS) l/2 , z max = Dsin(fo max ) ~ Db mm 
and z min = D sin(fo min ) ~ Db min - 

From eqs. (7) and (8) we can define a generalized Schmidt's 
(1968) estimator for the LF: 



where the sum is over all sources with luminosity within the 
interval AlogL around L. In eq. (8), for the vertical density 
profile we use a function of the form: 



p(z\R) = exp[-0.5(z/cr z ) 2 



(10) 



If Zmax is much larger than the vertical scale length a z and 
Zmin = -Zmax, the integral over dz in the right-hand side of 
eq. (8) gives {2i:) l l 2 a z . 

Bronfman et al. (2000) have investigated the scale height 
distribution of a population of 748 OB stars, for which (their 
Table 2) they tabulate values at different Galactocentric radii, 
for the North and South hemispheres, both separately and 
combined. We have fitted their quoted values for the com- 
bined distribution, after correcting for Rq = 8 kpc. The data 
result to be well represented by a function of the form: 



Zi/2 = 2.85# L69 pc 



(11) 



where z\n = {2\n2) l l 2 a z . This relation is also consistent with 
the range of values (39 to 52 pc) quoted in Paper II for HII 
regions in our Galaxy, therefore we adopt it. 

The choice of the radial density profile is more complex 
(see Fig. 12 of Paper II). In § 5 we have derived the distri- 
bution in Galactocentric distances of the sources in our sam- 
ples located in the fourth (Fig. 3) and first (Fig. 4) quadrants. 
Based on these distributions, we model the n(R) as: 



n(R) = exp[-((R-R pkl )/a rl ) 2 ]- 
+A exp [-((/? - Rpid)/ cr ri) 2 



(12) 



In the expressions above, R p ti, Rpk2, 0>i, Pr2 and A are free 
parameters to be determined fitting the observed distribution 
of sources as a function of the galactocentric distance R, given 
by: 

/■0max(/?) 

Af(R)AR = n(R)RAR d6 ■ (13) 

dzp(z,R) / dlogL(p(L) 

We note that the adoption of a double-Gaussian profile for 
the first quadrant is motivated by the fact that the line of sight 
intercepts the Carina and Scutum-Crux spiral arms. 

To keep the number of parameters within a manageable 
limit, we adopt, in eq. (13), the usual power-law representa- 
tion for the luminosity function <fi(L)dlogL oc L" 1 <f logL (i.e. 
a = 2), after having checked that the results are little affected 
by changes of its slope within the currently accepted range, 
1.5 < a < 2.5. The estimate of 4>{L) is thus based on an iter- 
ative approach: we use a first guess on the shape of <fi(L) to 
derive the parameters for n(R) that we then exploit to obtain a 
better estimate for the LF. More details on the calculation are 
provided in the Appendix. 

The minimum \ 2 values of the parameters for both sam- 
ples are listed in Table 2 and 3. In view of a derivation of a 
LF for different levels of completeness, we have derived the 
best-fit parameters for the radial distribution for each value of 
flux cut-off. For the reference cases of flux cut-offs at 5 and 
2 Jy for the CH87 and Lockman's samples, we obtain that the 
best fits correspond to a x 2 of 2.08 and 1.13, for a bin size 
AR = 1 kpc. The value of \ 2 is actually an overestimate and 
does not denote an unacceptably poor fit because the adopted 
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TABLE 2 

Best-fit parameters of the radial profile for the 4th 
quadrant 



Parameter 


5 Jy 


7Jy 


9 Jy 


Rpkl 


4.28 (kpc) 


4.26 (kpc) 


4.27 (kpc) 


Rpia 


9.54 (kpc) 


9.53 (kpc) 


9.53 (kpc) 


cr H 


1.44 (kpc) 


1.45 (kpc) 


1.45 (kpc) 


°72 


0.42 (kpc) 


0.42 (kpc) 


0.42 (kpc) 


A 


0.30 


0.29 


0.30 



TABLE 3 

Best- fit parameters of the radial profile for the 1st 
quadrant 



Parameter 


2Jy 


3 Jy 


4 Jy 


Rpkl 


5.17 (kpc) 


5.17 (kpc) 


5.16 (kpc) 


Rpk2 


3.70 (kpc) 


3.71 (kpc) 


3.69 (kpc) 


c r \ 


0.83 (kpc) 


0.82 (kpc) 


0.82 (kpc) 


a r2 


0.39 kpc 


0.39 (kpc) 


0.41 (kpc) 


A 


6.50 


6.51 


6.68 



TABLE 5 

Best- fit values for the LFs in the 1st quadrant for different 
levels of completeness (single power-law). 



Parameter 


2Jy 


3Jy 


4 Jy 


a 

x 2 
x 2 


1.85±0.11 
11.97 
1.33 


1.76±0.14 
23.02 
2.56 


1.88±0.12 
11.53 
1.28 


TABLE 6 

Best-fit values for the LFs in the 4th quadrant for different 

LEVELS OF COMPLETENESS (TWO-COMPONENT POWER-LAW). 


Parameter 


5 Jy 


7Jy 


9 Jy 


a 1 

logLknee 

logA 

x 2 
x 2 


2.02±0.11 
2.27±0.19 
23.45±0.11 

2.13 

13.43 

1.22 


2.03±0.15 
2.41±0.20 
23.54±0.08 

2.19 

9.88 

1.09 


1.74±0.22 
2.51 ±0.20 
23.47±0.12 

2.04 

8.69 

1.08 



TABLE 4 

Best-fit values for the LFs in the 4th quadrant for different 
levels of completeness (single power-law). 



Parameter 


5 Jy 


7Jy 


9 Jy 


a 


2.23±0.07 


2.22±0.09 


2.13±0.12 


2 

X 


19.49 


17.34 


27.81 


x 2 


1.49 


1.58 


2.78 



errors correspond only to the Poisson fluctuations of the num- 
ber of sources in each bin and do not include the contributions 
from the uncertainties in the estimates of distances R. Thus, 
although we are well aware that the axially symmetric distri- 
butions of eq. (12) are certainly an oversimplification, a more 
refined model, entailing more parameters, does not appear to 
be warranted in the present data situation. 

8. THE LF IN THE FOURTH AND FIRST QUADRANT 

The LFs for the CH87 and Lockman's samples derived by 
using eq. (9) and the radial profile n(R) determined in the 
previous section are shown in Fig. 5 and 6. The bin size is 
AlogL = 0.2. This choice makes it possible to perform di- 
rect comparisons with LFs obtained, with the same binning, 
by other authors (e.g. Smith & Kennicutt 1989 or KEH89). 

The errors in each bin are computed by applying the for- 
mulae by Gehrels (1986) (see § 3). In this case, the effective 
number of sources per bin, n e fr, is given by: 



»eff : 



(Si V m ) 



(14) 



A single power-law fit to the data retrieves the slopes pro- 
vided, for each quadrant and flux cut-off, in Tables 4 and 5. 
Quoted errors are retrieved from the fitting procedure and are 



TABLE 7 

Best- fit values for the LFs in the 1st quadrant for different 
levels of completeness (two-component power-law). 



Parameter 


2Jy 


3 Jy 


4Jy 


a' 


1.71±0.12 


1.68±0.14 


1.96±0.15 




1.98±0.33 


2.22±0.35 


1.93±0.37 


logLknee 


24.05±0.23 


23.94±0.45 


23.11±0.53 


logA 


1.75 


1.84 


1.88 


x 2 


15.47 


21.47 


14.24 


x 2 


2.21 


3.07 


2.03 



only indicative of the actual errors. We have also tried to fit 
the data with a two-component power-law of the form: 



4>(L)=A- 



(^/^knee) if L <C ^knee 

(L/Lknee) ^ if L > Lknee 



(15) 



From Table 6 and 7, we see that, for the CH87 sample, 
since A% 2 > 2.3 (see, e.g., Bevington & Robinson 1992), 
the addition of two parameters in the fit is justified and the 
two-component power-law model appears to be favored for 
all completeness levels, while for the Lockman (1989) sam- 
ple, the single power-law model always corresponds to the 
best x 2 value. 

To facilitate comparisons with published LFs, we convert 
our monochromatic LF (in units of ergs" 1 Hz -1 ) into Lyman 
continuum fluxes, Q, (in photons s" 1 ) and into Ha luminosi- 
ties, L(Ha), (in ergs" 1 ). To convert into Lyman continuum 
fluxes, we make use of eq. (6). 

To express the LF in terms of Ha luminosities, we use the 
relation (Osterbrock 1974): 



L(H(3)/hv Hfi a%{T e ) 



Q 



(16) 



which implies that the number of photons emitted by the neb- 
ula in H/3 is directly proportional to its Lyman continuum flux. 
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FIG. 5. — LFs in the 4th quadrant for different completeness levels. The solid line indicates the best-fit LF. Also shown (dashed line) is a two-component 
power-law fit. 



The hydrogen recombination coefficients and ap are tab- 
ulated, for various electron densities and temperatures, by Os- 
terbrock (1974) and Storey & Hummer (1995). By assum- 
ing an electron density n e = 10 2 cm" 3 , we derive the recom- 
bination coefficients at T e = 6200 K (CH87) and T e = 7400 K 
(Lockman 1989) by interpolation over the tabulated values. 
We adopt the value of 4.08 eV for the energy of a photon emit- 
ted in a transition between quantum states n = 4 and n = 2. The 
Ha luminosities are then obtained from L(Hf3) by interpolat- 
ing at 6200 K and 7400 K the Balmer decrements compiled 
by Osterbrock (1974) for different values of T e (Figure 7 and 
8). 

8.1. Effect of the choice of the radial density profile 

In this subsection, we discuss how the LF determined 
through application of eq. (9) depends on the adopted radial 
density profile. To illustrate the effect of n(R), we compare 
the LFs obtained for the CH87 sample by using three differ- 
ent functional forms, namely: the radial profile of eq. (12) 
with the best-fit values for the parameters in Table 2, an ex- 
ponential profile, a constant profile. The reference case for a 
flux cut-off at 5 Jy is considered. We chose the exponential 
radial profile derived by Kent, Dame & Fazio (1991) from the 
analysis of the Infrared Telescope (IRT) 2.4/im data: 



mean of the values tabulated by Bronfman et al. (2000) for 
the surface density of OB stars in annuli of increasing Galac- 
tocentric radius (see Table 2 of their paper): 



n(R) = 1 source kpc 



(18) 



n(R) = exp(-R/R d ) 



(17) 



with R d = 3 kpc. The expression above is used here only for 
illustrative purposes, therefore a discussion regarding the ex- 
act value of the radial scale length, Rd, is beyond the scope of 
the present paper. In the case of a constant profile, we use the 



Once again, the value of the constant is only indicative. Fig. 9 
shows the LFs resulting from using, alternatively, a double- 
gaussian [eq. (12)], an exponential [eq. (17)], and a uniform 
[eq. (18)] radial profile. Table 8 provides additional infor- 
mation, i.e. the best-fit values of the parameters (power-law 
index and break luminosity) for each curve. We have used a 
two-component power-law model given that, in the previous 
section, we have shown that it better represents the data for 
the CH87 sample. 

Clearly, the choice of the density profile affects the estimate 
of the luminosity function, although the derived values of the 
parameters are consistent with each other, within the errors. 

The different normalizations for the 3 cases reflect the dif- 
ferent effective volumes associated to each density profile. 
This is also clearly illustrated by eq. (7), where one can see 
that the effective volume acts as a normalization constant, 
while the shape of the LF is determined by the observed quan- 
tity N{L)d\ogL. 

8.2. Effect of the choice of the rotation curve 

In principle, the LF derived in § 8 could be dependent on the 
choice of the rotation curve model applied to the data to con- 
vert the measured line velocities into Galactocentric and solar 
distances. As discussed in § 5, our adopted rotation curve 
model is the one by FBS89. To illustrate the possible depen- 
dence of our results on this choice, we recomputed the LF for 
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35.4 
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Loq L(Ha) [erq*sec '1 
36.# 37.4 ^ 37,4 



36.4 
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7.4 37.b 38.4 
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. 4 h 

5 3 h 
* 2^ 
1 j- 

oL 

47 



49 

Log Q [sec" 



FIG. 7. — LF in the 4th quadrant expressed in units of Lyman continuum 
flux (lower scale) and Ha luminosity (upper scale). Results refer to a 5 Jy 
flux cut-off. 

the CH87 sample by using two alternative models, i.e. the 
rotation curve by Clemens (1985) and the one by Brand & 
Blitz (1993). The comparison between the LFs obtained by 
adopting these models with the one derived by using FBS89 
is shown in Fig. 10. Clearly, the three LFs almost overlap 
with each other, showing that our results are not sensitive to 
different choices for the rotation curve. 

8.3. Effect of the inclusion of sources with \Vlsr\ < 10 km/s 

In this section, we investigate how the removal of sources 
with | Vlsr I < 10 km/s affects the derived LFs. At this purpose, 
we recompute, for each sample and considered completeness 



r 



-1 Lj_ 

48.0 



48.5 



49.0 



49.5 
Log Q [sec"' 



50.0 



50.5 



51.0 



FIG. 8. — LF in the 1st quadrant expressed in units of Lyman continuum 
flux (lower scale) and H a luminosity (upper scale). Results refer to a 2 Jy 
flux cut-off. 

level, the effective volume, after examination of the radial dis- 
tribution of the sources. These are illustrated in Fig. 3 and 4. 
As noted in § 5, the removal of the low-velocities sources in 
the fourth quadrant causes a dip in the distribution for R ~ 8 
kpc. For this reason, instead of applying equation (12) to fit 
for n(R), we use a gaussian profile with inner cut-off, of the 
form: 



n 4th (R) = exp [-((R-R pk )/a r ) 2 ] x exp [-((R -R int )/R)\\9) 

For the Lockman (1989) sample, as in § 7, we adopt 
eq. (12), since the addition of the sources with \Vlsr\ < 10 
km/s leaves the radial distribution basically unchanged. For 
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TABLE 8 

Best-fit values parameters of the LFs for the 4th quadrant 
(5 Jy flux cut-off) obtained with different radial density 
profiles 



n(R) 


a'(L< L kncc ) 


P'{L > L kllcc ) 


logZ-knee 


eq. (12) 


2.19±0.20 


2.46±0.0.15 


23.45±0.11 


exponential 


1.99±0.19 


2.06±0.42 


23.47±0.09 


constant 


1.92±0.12 


2.98±0.63 


23.33±0.15 



NOTE. — Fit for a two component power-law LF [eq. I15H 



TABLE 10 

Best-fit values for the LFs in the 4th quadrant for different 
levels of completeness (single power-law). sources with 
\v lsr \ < 10km/s are also included. 



Parameter 


5 Jy 


7Jy 


9 Jy 


a 


2.09±0.06 


2.14±0.06 


2.17±0.08 


x 2 


26.96 


20.55 


22.40 


2 
Xr 


1.58 


1.21 


1.49 



35.4 35.9 
6 - 



Loq L(Hct) [eru*sec '] 
36.4 36. $ 37.4 37.9 



38.4 38.9 39.4 



2 - 



* n(R) from eq (12) 
o n(R) = exp(-R/R_d) 
A n(R) = const 



I* 



47 48 49 50 51 

Log Q [sec"'] 

FIG. 9.— LF for the CH87 sample (5 Jy flux cut-off) obtained by adopt- 
ing: a double-gaussian radial profile (eq. (12), stars), an exponential profile 
(diamonds), a constant profile (triangles). 



6 

i— , 5 



FBS rotation curve X 
Clemens rotation curve O 
BB rotation curve A 



22 23 24 

Log L [erg/(sec*Hz)] 



FIG. 10. — Comparison between LF obtained for the CH87 sample by 
adopting different rotation curves. A flux cut-off of 5 Jy is adopted. 



TABLE 9 

Best-fit parameters of the radial profile for the 4th 

QUADRANT. ALL SOURCES (INCLUDING THOSE FOR WHICH \Vls r \ < 
10KM/S) ARE CONSIDERED. 



Parameter 


5 Jy 


7Jy 


9 Jy 


R„k 


3.21 (kpc) 


2.82(kpc) 


3.13 (kpc) 


Rim 


1.94 (kpc) 


2.08 (kpc) 


1.99 (kpc) 


o> 


2.85 (kpc) 


3.24 (kpc) 


3.09 (kpc) 



illustrative purpose, we provide the best-fit parameters for the 
radial profile (Table 9), as well as the slopes (Table 10) of 



the LFs (Fig. 10), for the fourth quadrant only, given that the 
CH87 sample is also the on more impacted by the presence 
of sources with low measured velocities. From inspection of 
Table 10 and Fig. 10 and comparison with Table 4 and Fig. 5, 
we see that the slopes are consistent with the values retrieved 
by removing from the sample the low-velocity sources. In the 
case considered here, though, we find no evidence for a lumi- 
nosity knee. However, such a result should be interpreted with 
caution, due to the large uncertainty associated to the sources 
with \V LS r\ < lOkm/s. 

9. DISCUSSION 

In the following, we will refer, unless stated otherwise, to 
the reference cases of a flux cut-off at 5 Jy for the CH87 sam- 
ple, and at 2 Jy for the Lockman (1989) sample. In addition, 
we consider the case for which sources with \Vlsr\ < lOkm/s 
have been removed from the samples. 

In the fourth quadrant, the Ha luminosities of the CH87 
sample span ~ 3 orders of magnitude, from 10 36 to 10 39 
erg/sec or, in equivalent Lyman continuum luminosities, from 
1 47 6 to 1 50 6 photons s" 1 . The lower limit corresponds to HII 
regions ionized by a single star, while the upper value is in- 
dicative of the lack, in our Galaxy, of supergiant HII regions 
such as 30 Doradus in the Large Magellanic Clouds, which 
originate from very large OB associations. 

The LF derived from the CH87 sample presents marginal 
evidence for a change of slope at Lknee ~ 10 23 45 ergs" 1 Hz , 
corresponding to an Ha luminosity of ~ 10 37 75 erg/sec and 
to a Lyman continuum luminosity of ~ 10 4941 photonss _1 
(Fig. 7, upper and lower scales). A break in the LF of HII 
regions has been observed by many authors in external galax- 
ies. Among these, KEH were the first to conduct a systematic 
study of the HII region LF in 30 spiral and irregular galaxies. 
For 6 galaxies (of Hubble type Sab-Sb) of their sample, they 
find that the LF has an abrupt turnover at log L(Ha) = 38.7- 
39.0. Recently, Bradley et al. (2006) have reported a mean 
value for the break luminosity of log L(Ha) = 38.6 ±0.1, 
based on a sample of 56 spiral galaxies. Such a break has been 
interpreted in several ways. KEH suggest that the change in 
the slope corresponds to a transition between normal HII re- 
gions (ionized by a single star or by a small association) and 
the class of supergiants. Oey & Clarke (1997) seem to favor 
the hypothesis that the break might be caused by evolution- 
ary effects and a maximum number of stars per cluster, while 
Beckman et al. (2000) argue that it is due to the fact that, at lu- 
minosities above a given threshold, density bounding prevails 
on ionization bounding. 

As for the slope of the LF, our findings are in agree- 
ment with published estimates from extragalactic studies as 
well as with results obtained for our own Galaxy. For in- 
stance, KEH report, for Sbc-Sc galaxies, a power-law index 
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FIG. 1 1. — LFs in the 4th quadrant for different completeness levels. The solid line indicates the best-fit LF. All sources (including those with measured 
\ v lsr\ < lOkm/s) are considered. 



in the range 1.5 - 2.5 7 , while McKee & Williams and Smith 
& Kennicutt (1989) find, respectively, a = 1.99 ±0.25 and 
a = 2.3 ± 0.5 by fitting Lyman continuum luminosities greater 
than 10 49,5 photonss _1 and using the SMB sample. These val- 
ues are both consistent with the slope obtained for a single 
power-law fit, a = 2.23 ± 0.07, as well as with /?' = 2.27 ± 0. 19 
that we quote in § 7 for L > Lknee, in the case of a double 
power-law. 

In the first quadrant, the Ha luminosities (or the equivalent 
Lyman continuum luminosities) of the Lockman et al. (1989) 
span a smaller range of values than the CH87 data set, i.e. 
from 10 37 to 10 39 erg/sec or, in equivalent Lyman continuum 
luminosities, from 10 48 6 to 1 50 6 photons s" 1 . In this case, 
the LF is significantly flatter than in the fourth quadrant, sug- 
gesting a possible intrinsic difference between the HII regions 
population in the two quadrants, as also hinted by the discrep- 
ancy in the radial profiles (see § 4.). 

The LF obtained for the fourth and first quadrants were used 
to estimate the total ionizing luminosity of Galactic HII re- 
gions, given by: 



TABLE 11 

Total luminosities derived from the LFs in the 4th quadrant 



^tot = Vgp 



n(L)LdL 



(20) 





5 Jy 7 Jy 


9 Jy 


V G /.(kpc 3 ) 
LaalXHz (erg/sec Hz) 
L Ga l,Q (sec" 1 ) 
Laal,H a (erg/sec) 


8.98 8.85 
0.95 X 10 27 1.30X 10 27 
0.87X 10 53 1.19X 10 53 
0.19x 10 42 0.26x 10 42 


8.99 
0.63 x 10 27 
0.58 x 10 53 
0.13x 10 42 


TABLE 12 

Total luminosities derived from the LFs in the 1st quadrant 




2 Jy 3 Jy 


4 Jy 


V G /.(kpc 3 ) 
LaalSGHz (erg/sec Hz) 
L Ga i,Q (sec" 1 ) 
L Ga l,H a (erg/sec) 


10.49 10.77 
1.1 lx 10 27 0.82x 10 27 
0.94x 10 53 0.70x 10 53 
0.21 x 10 42 0.15x 10 42 


10.45 
1.31 x 10 27 
l.llx 10 53 
0.24X 10 42 



where n(L) = <fi(L)/(L In 10) and Vgp is the volume of the 
Galactic Plane occupied by classical HII regions, i.e.: 



y GP = 47T 



n(R)RdR / p(z)dz 



(21) 



7 Retrieved from single power-law fits. 



with n(R) given by eq. ( 1 2) and the best-fit parameters in Ta- 
bles 2 and 3, and p(z) as in eq. (10). In addition, R mdx = 15kpc, 
and, for the fourth quadrant, L m ; n = 10 218 erg/sec Hz, L max = 
10 25 erg/sec Hz, while, for the first quadrant, L min = 10 225 
erg/sec Hz, L max = 10 24 5 (see Fig. 5 and 6). The total lumi- 
nosities obtained by application of eq. (20) are shown in Ta- 
bles 1 1 and 12. Also provided are the values derived for flux 
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completeness levels above the reference cut-off. We note that 
our estimates are lower limits to the actual total luminosities, 
given that eq. (20) is integrated between a limited luminosity 
range. From the values quoted in Table 1 1 and 12, the average 
of the fourth and first quadrant is 0.89 ±0.23 x 10 53 sec _1 . 

The total HII regions Lyman continuum flux for the refer- 
ence cut-offs can be compared to the total ionizing flux of the 
Galaxy. Bennett et al. (1994) find Q = 3.5 x 10 53 sec _1 , based 
on COBE observations of the [Nil] 205 [im line, while McKee 
& Williams (1997) quote a value of Q = 2.56x lO^sec" 1 . 
Recently, Murray & Rahman (2009), from the analysis of 
WMAP data, have estimated Q = 3.2 x 10 53 sec _1 . Remark- 
ably, these numbers account for the combined contributions 
of individual HII regions and the diffuse warm ionized gas 
(WIM). In particular, HII regions are likely responsible for 
only a fraction (between 10 and 30%) of the total emission, 
as found, among other authors, in our analysis of the 5-GHz 
Galactic Plane free-free emission (Paladini, deZotti, Davies 
& Giard, 2005). Of the HII regions contribution, a significant 
amount is supposed to originate in sub-giant, or small HII re- 
gions, as noted, for instance, by McKee & Williams (1997). 
In particular, giant and sub-giant HII regions are likely to con- 
tribute 20% of the total ionizing luminosity of the Galaxy, 
while another ~ 20% is supposedly due to small HII regions 
(see their Table 3). McKee & Williams (1997) find a total 
ionizing luminosity from giant and super-giant HII regions of 
0.57 x 10 53 sec _1 which, scaled by a factor accounting for the 
difference in the adopted Rq value (8.5 in their case), becomes 
equal to 0.50 x 10 53 sec . Taking into account the contribu- 
tion from small HII regions, this number needs to be doubled, 
to roughly 1 x 10 53 sec _1 , i.e. one third of the total Lyman 
continuum flux of the Galaxy, consistent with our results. 



10. CONCLUSIONS 

We have investigated the LF of Galactic HII regions in the 
fourth and first quadrants, using the samples from the CH87 
recombination lines survey for the fourth quadrant, and from 
the Lockman (1989) recombination lines survey for the first 
quadrant. We have shown that a good determination of the 
highly inhomogeneous spatial distribution of HII regions is 
essential to obtain an accurate estimate of the LF: commonly 
used but inaccurate density profiles affect both the shape and 
the normalization of the LF. 

The LF in the fourth quadrant is well-represented by a 
single power-law with spectral index a = 2.23±0.07, al- 
though we find that \ 2 goodness-of-fit analysis slightly fa- 
vors a two-component power law with a break luminosity at 
log[L(//o;)/erg/s] = 37.75. The LF in the first quadrant ap- 
pears to be flatter, and the best-fit is obtained in this case for a 
single power-law with slope a = 1.85±0.11. 

The inferred total contribution of individual HII regions to 
the global ionizing luminosity of the Galaxy is found to be in 
agreement with the estimates by McKee & Williams (1997) 
and roughly equal to one third. 
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APPENDIX 



From eq. R = (flg + D 2 -2D/? cos0 1/2 we obtain: 



D = R cos(/) ± (R 2 - R 2 Q sin 2 (0) 



1/2 



(1) 
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If we call 8 the angle between the direction from the Galactic Center (GC) to the Sun and the direction from the GC to the source 
(i.e. between Ro and R), the sine theorem gives: 

sin(0) = ^sin(/). (2) 



Combining these equations we get: 



sin(6>) = 



cos(/)± 1- 



1/2" 



sin 2 (/) 



sin(Z) . 



The number of sources detected at distance R, within AR, from the GC can be written as: 



t*9m' dX (R) r-Zn 

Af(R)AR = n(R)RAR / dO \ 



(R,9) f logL 

max 

dzp(z,R) / dlogL cj)(L) 



(3) 



(4) 



>0mm(R) Jzmm(R,0) J logL mi „(R,6>) 

where z mdX Mn(R, 0) = b mdx . min D(R, 9) and L min (R, 9) = 4irD 2 (R, 9)S lim , with D(R, 9) given by: 

D = (Rl + R 2 -2RR cos9) l/2 . 

If we have a sample selected in the longitude range / m ; n < / < l mdx and R<Rq sin l m . dx (i.e. R is smaller than the minimum distance 
between the GC and the line at I = l mnx , assuming I < 90°), we have: 



(5) 



sin 9 mdxMn (R) = { -^cos l mia ± 



1 - I ) siir /min 



1/2- 



• sin L 



(6) 



where 9 mdx obviously corresponds to the + sign. Note that 9 mdx > ir/2, so that 9 m . dx = tt - arcsin(sin 8 m . dx ), while 9 n 
arcsin(sin # m ; n ). 

If R > Rq sin l mdx , we do not have contributions to Af(R) from the range 9\ <9 < 92, given by: 



sin6» 2 .i(7?) = { yCos/,„ : , x : 



sin L 



1/2' 



sin /„ 



(7) 



In practice, we split the integral over 9 in 2 parts, from 6Vn to 9\ and from 9^ to 9 m . dx , where 6*2,1 are given by the eq. (A7) if 
R > Rosinl max or by 9\=9i = Tt/2-l mdx otherwise. 



